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ON THE THEORY OF THE RECTILINEAR OSCILLATOR. 1 

By Edwin Bidwell Wilson. 
Introduction. 

Planck's rectilinear oscillator is of great importance in theoretical 
physics. To grant this it is not necessary to assert a belief in the 
physical existence of the oscillator, much less to assert Planck's 
belief in such existence. The Carnot engine may not exist; but the 
concept has proved so fruitful as to be fundamental in the general 
field of thermodynamics. The (perfect) semipermeable membrane 
may be a figment of the imagination; but, in special fields in thermo- 
dynamics, reasoning based on the membrane is of crucial impor- 
tance. And in that branch of thermodynamics which deals with 
radiant energy, the oscillator has held the position of cornerstone. 
Moreover, as in the case of the Carnot engine and semipermeable 
membrane, there does exist in physics a (more or less imperfect) 
prototype of the oscillator — the radiator or resonator of Hertz and 
his followers in radio telegraphy. 

In addition to the thermodynamic importance there is an interest 
attaching to the oscillator as the simplest model of the atom. Hypo- 
thetical models of the atom vary greatly; some are founded on 
statical configurations, others on dynamic equilibrium, some call for a 
number of corpuscles per atom comparable with the atomic weight, 
others for a number far greater than this, some represent the positive 
charge as widely spread out with the negative corpuscles imbedded in 
it, others concentrate the positive charge at a point. It is not neces- 
sary to take the oscillator very seriously as a model of an atom; all 
that is useful is to observe that by its means energy is absorbed and 
emitted and that thus it performs the functions of matter in establish- 
ing a temperature equilibrium in a closed field of radiant energy. 

The object of this article is to make some comments on the equation 
which defines the motion of the oscillator. These comments are 
obvious at first sight to any one who examines the differential equa- 
tions used in the theory of the oscillator, but as I cannot find them in 
the literature, I venture to print them at this late day. 

1. The differential equations. Let us consider first some of the 



1 Read to the Academy at its Stated Meeting, March 11, 1914. Received 
June 16, 1914. 
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prime facts or hypotheses underlying the analysis. These supposi- 
tions are as follows: 

1°. The oscillator executes in the main a simple harmonic motion. 
This leads to the statement of the (approximate) equation 

% + &x = 0, (1) 

where z is the coordinate defining the state of the system. 

The basic reason for this assumption is that it is desirable for the 
oscillator in its function of emitter of energy to be as simple as possible, 
that is, to emit monochromatic radiation of a definite frequency; 
for our customary analysis resolves heterogeneous radiation into 
monochromatic elements. 

(If we wish to assume that the oscillator consists of a widely spread 
uniform positive charge at rest and of a concentrated negative charge 
vibrating to and fro through the center of the positive charge, we may 
obtain a model of the oscillator which satisfies the condition that the 
force on the moving particle shall be proportional to the displacement ; 
for if a particle moves within a uniform sphere and if the elemental 
attraction follows the Newtonian law, the resultant force is proportional 
to the distance of the particle from the center of the sphere.) 

2°. The motion of the oscillator is not affected by frictional re- 
sistance. This leads to (or corresponds with) the omission of a possi- 
ble corrective term proportional to the velocity dx/dt in (1). 

Perhaps the best reason for omitting frictional terms is that there 
seems to be no necessity for complicating the mathematical equations 
or physical hypotheses by admitting friction. In the Hertz radiator 
there is frictional dissipation of energy due to the resistance of the 
circuit; in the idealised oscillator the resistance would naturally be 
omitted. 

3°. The motion is damped by the emission of energy. If the form 
of (1) appropriate to the principle of energy be used, that is, the 
integrated form, the undamped motion would be characterized by the 
equation 

fdx\ 2 



§W=C or ± 



, , +k 2 x 2 
dt 



0, 



and the damping would then be inserted by a change to the form 

dx\ 2 
dt) Jo dt 

where F is the rate of radiation. 



dt < + kV 



- 2F, (2) 
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The modification of the equation is thus based directly on the 
principle of energy. The form in which F is thrown depends some- 
what on the point of view. 

If the oscillator is considered as a Hertz dipol and the ordinary 
solution for the electromagnetic field set up by this radiator at large 
distances from its center is obtained, we find by the application of 
Poynting's theorem that a certain amount of energy is radiated off 
per period and, in particular, that the amount is proportional to the 
square of the acceleration. Then the equation, written for a com- 
plete integrated period T is 



f 

J to 



h+T 



+ k 2 x 2 



2m 



<PxV 

dty 



dt= 0. 



An integration by parts transforms the third term. Whence 






+ T 



+ kW 



-2m 



dxd 3 x 



dt d? 



dt 



2» 



«o+ T 



= 0. 



(3) 



fdx\ 
\dtj 
dx d?x 

jtd¥_ 

Now if to be an epoch for which either the velocity or the accelera- 
tion vanishes, the last term drops out. Under these conditions 

a + t { §_ r fdx\ 

J to \dt[\dtj 

The conclusion is then drawn that 



+ k 2 x 2 



„ tit// (1 *C i -j 

~ 2lx dtdF Ut: 



o. 



d 

dt 



[ 



f) +kV 



2a 



dx d s x 
dtd¥ 



= 0. 



(4) 



Exceptions can be taken to this conclusion whether the argument is as 
above or as given in slightly different form by Planck. 2 It by no means 
follows that, because the integral of a function over a complete period 
from tn-\-nT to to-\-(n-\-l)T vanishes, the integrand must therefore 
be identically null. 

What we may infer is that 



+ kV 



■*■££-» 



where a(t +(n+l)T) — a(t +nT) = 0. There is nothing at all to 



2 Planck, Vorlesungen iiber die Theorie der Warmestrahlung, first edition, 
110 (1906). Abraham, Theorie der Elektrizitat, second edition, 70 (1908). 
Rudenberg, Annalen der Physik, 25, 346-466 (1908). (In his second edition, 
Planck introduces the quantum hypothesis early and omits this derivation of 
the equation of the oscillator.) 
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suggest that the function a(t) has a negligible effect upon the motion. 
It therefore follows that (4) or the equation 

^ + Fa--M~=0, (5) 

obtained therefrom by discarding the factor dx/dt, is at best only 
approximate. Concerning the degree of approximation we may 
perhaps be willing to admit that the error is always extremely small ; 
but we should be very chary about admitting that the error is neces- 
sarily small relative to the terms retained. 

That (5) is only an approximate result is clearly brought out by 
Planck in his derivation, but the clearness is somewhat obscured by 
a footnote referring to Abraham's derivation in a manner which 
implies, though it does not state, that the latter proof is exact. Abra- 
ham gives in the first instance essentially the proof which we have 
given above. He goes on, however, to show by a discussion of 
the physical dimensions of the terms that the equation (5) is the 
only possible one. What he really proves is that under certain hy- 
potheses (5) is the only possible linear equation of the form 

tf< 2 + i °< dt 

The physical dimensions are equally well satisfied by the non-linear 
equation 



dx 
~dt 



g + H* 



+ M(g)=<>, (0) 



which follows naturally and immediately from (2), that is, from the 
principle of energy, provided it is assumed that the rate of radiation 
is proportional at each instant to the square of the acceleration, instead 
of merely proportional thereto when integrated over a complete 
period, and which would arise from the equation in integrated form, 
provided one were willing to apply to it directly the reasoning which 
is applied after an integration by parts. 

Lorentz has also derived (5) by a method which throws emphasis 
on the approximate nature of the result. 3 That he regards the equa- 
tion with some suspicion may be inferred from his statement of a 
later date 4 : Eine Gleichung "die iibrigens den Mathematikern noch 
manche Frage darbote." 

3 See, for example, his Theory of Electrons, 251, (1909). 

4 Physikalische Zeitschrift, 11, 1250 (1910). It was this remark which so 
strengthened my own suspicions as to lead me to take up the question some 
time ago. 



WILSON. — RECTILINEAE OSCILLATOR THEORY. Ill 

The calculation of the instantaneous rate of radiation of energy 
from a moving electron, as contrasted with the integrated rate hitherto 
used, has been carried out by a number of authors in different ways. 
The rate is almost universally found proportional to the square of the 
acceleration when the square of the velocity is negligible relative to 
the square of the velocity of light. 5 The mathematical statement 
of this result is contained in (6), and this equation therefore might 
more properly be taken as representative of the motion of the oscil- 
lator than the equation (5), — provided, of course, that we are willing 
to use the electron model of the oscillator. 

A comparison of (5) and (6) will throw some light on the legitimacy 
of the approximations used in deriving (5). The terms to be com- 
pared are 

dx cPx , /dV\ 2 

— M ^t ~77J and » l — I 



dt df ^\dfij 

On the supposition that the oscillator executes in the main a simple 
harmonic motion, the signs of the two terms are the same. But the 
phases of the magnitudes are very different; for the first term is a 
maximum as the particle passes through the position of equilibrium 
and vanishes in the positions of extreme elongation, whereas the 
second term is a maximum in the extreme positions and vanishes in the 
position of equilibrium. The average magnitudes of the terms are 
the same, but it is difficult to see how two terms of equal average 
magnitude could be more different than these two in their effect upon 
the equation of motion. 

We shall now discuss the two forms (5) and (6), the first rather 
briefly, because it has been tolerably well treated, the second more 
in detail, because it seems to have been neglected. 

2. Discussion of (5). As we have to deal with cubic equations 
and their (approximate) solutions, it is necessary to determine the 
order of magnitude of the coefficients and it is desirable to put the 
equations into their simplest form. If m is the mass of the vibrating 
particle, e the charge in electrostatic units, c the velocity of light, and 
if we adopt the customary coefficient for the radiation term, we have 



5 For two different lines of deduction see J. J. Thomson, Electricity and 
Matter, chap. 3 (1904), and Wilson and Lewis, The Space-time Manifold of 
Relativity, These Proceedings, 48, 389-507 (1912), with especial reference to 
p. 480. A. Macdonald in his Electric Waves, 75 (1902), finds a steady rate of 
radiation from the oscillator. 
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or 
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+ k 2 x 
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cPxV 
dt 2 ) 
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As numerical values we may assume 
e = 4.7 10- 10 , m = 9 10" 

Hence 

2 e 2 



c = 3 10 10 . 



3 mc 3 
Substitute t = ar and choose ct 2 fc 2 

d 2 z 



= 6 10- 2 



1. Then 



and 



dx 
dr 



dr 2 



d?x 



+ a;-6)kl0- 24 f^=0 



+ x 



+ 6*l^(g) 1 : 



0. 



Now fc = 2irc/X, where X is the wave length. The observed range of X 
in the spectrum is from 10~ 5 to 3 10 -2 with X about 5.5 10~ 5 in the 
center of the visible spectrum. The range for the coefficient of the 
last term in the equations is therefore from about 10 -7 in the ultra- 
violet to about 4 10 _u in the infra-red with 2 10 -8 for the value near 
the middle of the range of visibility. 

The numerical forms of the equations are therefore 



d 2 x , d?x 

d? + X - K d? = ° 



and 



with 



dx 


d 2 x . 


+■(£;-« 


dr 


i? + 'j 



(5') 



(6') 



10- 7 > k > 4 10-" 



or 



2 10- 



the latter being taken as a sort of standard value. 

(It is not necessary to make the special assumptions involved in 
the electron model of the oscillator if we desire merely to determine 
the order of magnitude of k. If (5') be replaced by 



d 2 x dx 

-T-. + k- f +x= 0, 

flr 2 dr 



(7) 
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which is nearly equivalent to it, since x and cPx/df are nearly equal, 
the amplitude of the motion is seen to suffer a reduction in the ratio 

e- inm : 1 

during n complete oscillations. We may now use experiments on the 
distinctness of interference fringes for different path-lengths to 
obtain an estimate for k, and we find that 2 10~ 8 is an entirely 
reasonable value. See Abraham, loc. cit.) 

The complete solution of (5') will be of the form 

x=& e aT + c-" T (C 2 cos7T 4- C 3 sinYr), 

where a,—(H±yi are the roots of the cubic 

r 2 + r — k r 3 = 0. 

The root a is positive and of great magnitude, approximately 1/k, 
say, \ 10 8 . As Planck says, this root has no physical significance. 6 
This means that for every determination of the constants Ci, C, 2 C 3 , 
the value of C\ must vanish and the solution of the equation reduce 
to 

x = c~ &T (C 2 cos7r + C 3 shryr). 

It is certainly unfortunate that we are forced to discard so arbi- 
trarily one root of an equation which we have labored so carefully to 
establish. Moreover, as the solution of the equation now depends 
only on two constants, the elimination of these constants by differ- 
entiation will give an equation of the second order which accounts 
satisfactorily for the phenomenon as far as its mathematical side is 
concerned. This equation, which is almost identical with (7), con- 
tains a friction term, it is true, and is unsatisfactory physically be- 
cause of this fact. But there is no real reason for considering it as 
less satisfactory than (5') which has a physically meaningless root and 
which moreover, while ostensibly allowing for radiation on a physical 
basis, actually requires the rate of radiation to be maximum where 
there is a general agreement that it should be minimum, and minimum 
where it should be maximum. 

3. Discussion of (6). We turn next to the equation (6') and, 
for the sake of brevity, write it as 

K.f+vf+vx = 0, (8) 

/ and v designating respectively acceleration and velocity when the 

6 If we make clear the assumption of the quasi-stationary state as a prelimi- 
nary condition in our approximations which lead to the equation, the large 
root is sure to be inadmissible by virtue of our assumption. 
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unit of time is the fraction X/(2xc) of the second. This equation 
differs very minutely from (5') or (7), but there is no reason to believe 
that the motion defined by the equation should resemble very closely 
the motion defined by those. As the integration of the non-linear 
equation of the second order is not simple, we shall first examine the 
equation to find certain elementary qualities of the solution. 

A. When the velocity vanishes, the acceleration also vanishes. 

B. When the acceleration vanishes, either the velocity or the 
displacement vanishes. 

C. In the neutral position x — 0, either the acceleration vanishes 
or the acceleration and velocity satisfy the relation k/ + v = 0. 

D. Arbitrarily assumed (initial) values of the acceleration and 
position determine the velocity, and similarly values of the accelera- 
tion and velocity determine the position, but assumed values of the 
velocity and position allow the acceleration either of two values 

-_ v 1 , 

/ — — .7- ± K~ *V-—4:KVX- 

Ik Ik 

E. Values of v and x which make « 2 — 4kvx negative are not ad- 
missible. Hence if x is positive, v must either exceed 4xx or be nega- 
tive; and if x is negative, v must be positive or less than Akx (alge- 
braically) . 

F. The case of x positive and v greater than 4/ca; leads to a point 
d'arret. For under these conditions / is negative, no matter which 
sign is taken before the radical, and the velocity must be decreasing as 
x increases; when v becomes equal to 4/ca;, the acceleration becomes 
imaginary, and so does the motion. The same arises under the case 
of x negative and v less than 4kx. 

G. The case x positive and v negative also leads to a point d'arret. 
For under these conditions if the negative sign is taken with the 
radical, then / is negative and the magnitude of v is increasing as x 
diminishes toward zero, the particle will reach the origin with a definite 
velocity and pass through to the negative values of x where it falls 
under the case F; and if the positive sign is taken with the radical, 
so that / is positive, the particle is slowing down, and we have to 
consider two possibilities, first where the particle reaches the origin 
before its velocity vanishes and where the argument just given will 
then again be applicable, second where the particle comes to rest 
before reaching the origin only to start on again with a negative / and 
fall under the first supposition of this case G. A similar discussion 
may be given for the hypothesis x negative and v positive. (There 
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is a special case in which the particle just comes to rest at the origin, 
remaining there in a sort of unstable equilibrium.) 

The consequences of this discussion of the second type of equation — 
the equation (6) founded on the radiation formula generally adopted — • 
may be formulated in the unsatisfactory statement that, except for 
two special possibilities, the motion runs into a point d'arret where at 
a certain time and in a certain position, with a definite velocity and 
acceleration neither of which vanishes, the motion suddenly becomes 
imaginary. This is very disconcerting; it is at direct variance with 
our accepted notions of what happens in mechanical or electrodynamic 
systems. 

It is also a consequence that the oscillator governed by equation (6) 
cannot oscillate. 

Unfamiliar as this result may sound to those who have a feeling 
that a slight change in the differential equation introduces only a slight 
change in the motion, we must nevertheless admit the result and de- 
velop the contrary feeling that a slight change may entirely alter the 
type of motion. This fact has been brought out by Borel in a paper 
just printed 7 wherein it is shown that the equation 1? + x 2 = a 2 + H 
no matter how small the constant X, does not define an oscillatory 
motion, nor does the motion approach the oscillatory type as X ap- 
proaches zero. Any one might admit that, in the course of a long 
time, \t would become important and the type of motion be changed; 
but this is not the true explanation, for, no matter how small the 
constant X, the motion departs from the simple harmonic type from 
the time t = x/2, that is, after executing a quarter oscillation. 

4. The integrals of (6). We shall now turn to the integration 
of the equation 

dx (d 2 x , \ , f(Px\ 2 „ 
The first thing we shall treat is motion from a position of rest. In 



7 Borel, Annali di Matematica (3), 21, 225, (1913). The equation of 
Borel corresponds to a constant radiation of energy. The solutions of the 
equation do not exhibit the point d'arret, but are continuous. It has seemed 
to me, notwithstanding the publication of Borel's paper, and perhaps even on 
account of it, that the printing of my own results would be of interest. The 
comments on Borel's paper, which he prints as a footnote to a brief account of 
his results in his new book, Introduction Geometrique a Quelques Theories 
Physiques, 107 (1914) appear less applicable to the results of our present 
work. 
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this case v and/ are both zero. Let us therefore assume a solution for 
a; as a series 

x = o(l + ar 3 + j3t 4 + T t 6 + . . .) 

in r. The equations for the determination of the coefficients are 

36a 2 K + 3a = 0, 14Aa$K + 18a 2 + 4)3 = 0, 

(144(3 2 + 240a7)* + 60a/3 + 5? = 0, . . . . 



Hence 



or 



a = 0, = 0, 7 = 0,.... 



1 „ 1 11 



7 = 



12k' k 64k 2 ' ' 3840k 3 ' "" 

Thus one solution for x is x = a, the particle remains at rest in the 
initial position, and the other solution is 

/ T 3 T 4 11t 6 \ 

:E=O ( v 1 ~I^ + 64k 2 -38407 + ----> ' (9) 

the particle moves with an acceleration toward the origin, as was 
foreseen under G above. 

If this solution in series is valid, that is, if the series converges, it is 
at any rate of no value for purposes of calculation unless t is exceed- 
ingly small, that is, for intervals of time which are very small com- 
pared with the period of monochromatic light (t = 2ir). The value of 
the series is first to show how the motion starts and second to show 
that there is a possibility of motion, that is, that there is a solution of 
the equation other than the solution x = a with the given initial 
conditions. In discussing this question we shall follow the method 
given by Picard. 8 

Let the equation (6') of the second order be written as two simul- 
taneous equations of the first order, namely, 



-($h 


fdv 

'Kir 


+ z) = 0, 


dx 
dr 


-v= 


(10) 


These may be expressed 


in solved form as 










dv 
dr 


9 ± 1 
2k 2k 








dr 


V«2. 


— 4lKVX. 




The fundamental theorem on 


the existence of integrals 


fails for those 



8 E. Picard, Traite d' Analyse, 3, chap. 3. 
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values of the variables for which two values of the derivatives become 
equal. In this case the failure is when 



This equation factors, and there are therefore two surfaces 

v = 0, v = <Lkx 

in the space (v, x, t) which are singular in the sense that the funda- 
mental existence theorem fails. 

Of these surfaces one, namely, v = 0, satisfies the differential equa- 
tion (6') and is therefore singular also in the sense that it is the 
envelope of particular solutions. It is possible to put. the particular 
solutions in evidence. For let x be positive so that v <Q. Intro- 
duce a new variable v = — z 2 and let the sign of z be so chosen that 
z increases with t. We have then the equation 



dz 



:-r 



dr 2/c ' 

In cancelling the variable z the double sign goes out. For since / is 
decreasing with r (see G above) and since 

. dv „ dz 

ar dr 

it follows that dz/dt is positive. We therefore have to integrate the 
equations 

dz 2,1 ; dx 

dr = lK + ±K^ KX + Z ' dr=~ Z ' 

with the initial conditions 2 = 0, x = a < 0, to obtain the particular 
solutions which are tangent to the singular solution v = 0. For 
these equations the existence theorems are applicable. As a solution 
exists, it must be (9) which we found by undetermined coefficients. 

The other singular solution v = 4kx is not singular in the sense that 
it is enveloped by particular solutions but in the sense that it is a 
cusp locus for particular solutions (see Picard, loc. cit.). Without 
following Picard through the individual steps of the work we may 
write down the result. If v and x are two values satisfying the 
relation v = 4kx at r = t , the solution may be expanded into the 
form 

x = x + a (to - r) + (r - t)* + 7 (to - r) 2 + 5 (r - t)» + . . . 
v = Vo+ a' (to - r) + 0' (ro - t)» + y' (r - r) 2 + 5' (r - t)» + . . . 
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Now let us take the case a?o>0, /o = — Vo/2k. The constant /3 is zero, 
and 

x = x — v (r — t) - ^ (t - t) 2 + 5 (t — r) 1 + . . . . . , 
» = «o + £ (to - t) + - 5 (t - t) s + . . . 

The value of 5, which does not vanish, may readily be found, if desired. 

This solution brings out the fact very clearly that the motion runs 
into a point d'arret. The curve, to be sure, has a cusp, but the expan- 
sion is in powers of V To — t, which means that the time cannot exceed 
to. If we could imagine that time satisfies the same possibilities of 
reversal of direction as space, we could interpret this cuspidal motion, 
but as time can flow only in one direction, the motion must be inter- 
preted as having a point d'arret. This does not mean that the two 
parts of the geometric locus cannot be interpreted, but that they must 
be representative of two different motions terminating in the same 
point d'arret. 

We may indeed trace the motion back from the point d'arret along 
the two branches. We have 



f= -- 
•' 2k 



2k 



V« 2 — 4kto. 



The acceleration is negative for both signs before the radical; x has 
increased up to x at t = to and v has decreased to % In following 
the motion backward through time we find that x decreases and » 



X 




/A-) 






1 r 


-a 


(+)/ 





Figure 1. 

increases. For the motion with the positive sign before the radical 
the retardation is less than for the motion with the negative sign. 
Hence in going backward, the former motion gives larger values of x 
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and smaller values of v than the latter motion for the same value of 
tq — t. When the (+) motion passed the origin, the acceleration was 
zero; when the (— ) motion passed the origin, the acceleration was 
negative and twice as large as at the cusp. The (+) motion may 
be traced back to a point where the particle comes to rest for a certain 
negative value of x; this corresponds to the solution (9), and the 
motion may be traced still further back, the velocity and retardation 
increasing all the while. The (— ) motion may likewise be traced back 
through the origin x = 0, the velocity and retardation increasing. 
The curve which shows roughly the relation between x and r is 
given in Figure 1. If we call the curve which corresponds to the posi- 
tive sign x = o$(+)(t) and the curve which corresponds to the negative 
sign x = o$(_j(t), we may at once write the complete solution of the 
differential equation (6'), namely, 

x = &*, + ,(r - ft), x = Ci$,_,(r - ft), 

with the singular solutions x = C. The first solution could have 
been written at once from (9) ; the particle is at rest at x = d when 

T = "ft. 

The way in which the constants of integration enter into the solu- 
tion (which could be foreseen from the differential equation itself 
inasmuch as the equation admits the two-parameter group x' = x + a, 
t' — t + a) shows that the motion, though not periodic, has definite 
intervals of time independent of the constants of integration associated 
with it — the time from rest to the origin and the times (two possi- 
bilities) from the origin to the point d'arre't. 9 

5. Numerical Calculations. To ascertain these periods of time 
we must integrate the equation (6')- Let 
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dx 
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= — xu, 
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du\ 2 
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— u(u 2 - 


dT ) 
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= M 2 - 


u 1 

"2k ± 2k 










Vw 2 + 4km. 




Suppose x = 


■■ a <0, 


» = 


' / = 


are the initia 


1 conditions. 


Then u = 0. 


Moreover, as 


/ becomes positive, du/dr 


— M 2 


becomes 



1 These latter would not appear in the case Borel treated. 
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positive. Hence the positive sign must be taken with the radical. 
When x — 0, u has become positively infinite. The time of fall to 
the origin is therefore 



Introduce the substitution 



y 2 = 



W . J. , : . 



4kz/ 2 



u + 4k' 



1-2/ 2 



(13) 



Then 



Jo ■ 



4x dy 



&V + (l - 1/ 2 ) (i - y)' 

The integration may be performed by partial fractions. It is 
necessary to factor the denominator approximately. One root is 
nearly equal to —1, the other two are imaginary. The real root may 
be developed in powers of k 2 . The first approximation gives 

y = - 1 + 2k 2 . 
The denominator is factorable as 

(1 + 8k 2 ) (y + 1 - 2k 2 ) [y* - 2(1 - 5k% + 1 - 6k 2 ]. 
The reduction of the integrand to partial fractions gives 

1 y - (1 - 5k 2 ) - (2 - 7k 2 ) 



K(1 ~ 2K2) b/+l 
For the integration 



t = k (1 - 2k 2 ) 



i lo g; 



2k 2 y 2 - 2 (1 - 5k 2 ) y + 1 
(y + 1 - 2k 2 ) 2 



6k 2 



2 l " 6 ^ _ 2 (i _ 5 K 2) j/ + l - 6k 2 



, 2 - 7k 2 , - x w - 1 + 5/c 2 
+ ~^ tan 2k 



Substituting the limits and 1, we have 
t = k (1 - 2k 2 ) 



logl^^U 2 - 7 * 2 



1 -2k 2 



The k 2 is now everywhere negligible, and 



2k 



OK 



tan ~ + tan 



5k 2 \ 



2k ; 



r= 2 + K,0g K + 2- 



WILSON. — RECTILINEAR OSCILLATOR THEORY. 121 

As k is in the neighborhood of 2 10~ 8 , the difference between r and 7r/2 
is about 4 10 -7 . 

We desire next to follow the motion through the origin and on to 
the point d'arrSt. 10 Construct the Riemann surface for the integrand 
in (12). The surface has branch points at u — and u = — 4/c, and 
a junction line between them. We have let u trace the real axis from 
to oo ; now let u turn through a large circuit in the upper half -plane 




and come down to the other end of the axis of reals. No time elapses 
during this process, for approximately dr = du/u 2 . The variable x 
merely swings around the origin on an infinitesimal semi-circumference 
and is ready to start off on positive values as u comes in along the 
negative real axis. That the branch point u = — 4/c corresponds to 
the point d'arret for the motion is readily surmised, and may be 
proved from the defining equations (11). 

The time from the origin to the point d'arrSt is therefore 

,,_ 4k du 

T '= , « 

J— oo W — 



1 



2k 2k Vm2 + W 
But with the substitution (13) the integrand becomes again the same : 

r= r 4k dy 

Ji 8kY + (i - f) (i - y y 

Hence 

r-.o-*>[-*4*+4p:g-.i 1 f" 

10 The simplest method seems to be one dependent on the theory of functions 
of a complex variable; for it enables us easily to pass around the singularity 
at the origin. 
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Again k 2 is negligible and 

„, x 15 

If we had desired merely the total time T + T' = ir — 2k for the 
swing from rest to point d'arret, we could have obtained the result 
more directly by the theory of residues. The integral 

du 

T -4- 7" = / "2 u u I 4k> 

T+T Jo m2 -2,+ 2k\ 1 + w 

taken around the path in Figure 2 which passes through infinity, is 
equivalent to the sum of two integrals 

du „ du 

(14) 

U ±K 6K * ■ u 



~ du du 

J - 4k U 2k^2k\ 1 + 7, JoU 2k + 2k\ 1 + 77 



where the first is taken along the direct path on the upper side of the 
junction line and where the second is taken around that pole of the 
integrand which lies in the upper half of the upper sheet of the Rie- 
mann surface. As the calculation of the position of this pole is neces- 
sary for the work of the next section we shall now check the above 
value of T + T'. 
The denominator of the integrand vanishes when 

This gives one zero at the branch point u — 0; the integral to u = 0, 
however, is finite as the root is simple and the branch point has two 
sheets. The other roots are determined by the cubic 

— KU 3 + M 2 + 1 = = — K{U — ~ K — K)(u 2 + KU + 1). 

To terms in k the approximation gives the three roots 



1 . . K 

U= \- K, W==fcT — 

K 



9 



The first lies in the under sheet and does not concern us. The root 
u — i — \k lies in the upper half of the upper sheet, and the residue 
of the corresponding pole has to be found. The expansion 

2k 2k \ u u 2 J 
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of the denominator is legimate in the neighborhood of u = i — §k. 
The denominator therefore becomes 

— l+^Kjlu+l—- 

2 M 2mm 



1 1 t, -\ 



u< 


+ 1- 


K ' 
U 


The residue 


is 






1 


- 


u-\- i 


1 

~2 K 


K% 
U 



. 1 

u= 1 — 



2i (1 + vc) 2i 



2 



The second integral in (14) is therefore ir (1 —it), which is 2« times 
the residue. 
To calculate the first integral in (14) we may rationalize: 



-f" 

J— 4 K 



2k 2k\ 1 + „ 

du. 



•I 

K K 



. M° M 
M 4 



The first two terms of the denominator and the first of the numerator 
are of the order k 2 , or higher between — 4k and and may be neglected. 
The result is 



~X 4 ,2 + 2\I 1 + ^ M: 



-it 



u + V M (w + 4k) + 4k log ( V M + 4k + 5m) 
- 2k - 2k log (± i), 



o 

— 4« 



where the indeterminate sign and the indeterminate value of the 
logarithm must be properly chosen. This need not trouble us as the 
term is imaginary, and its only function is to cancel the imaginary 
part of x(l — in). The total value of T + T' is thus seen to be x — 2k 
as previously determined. 

We have next to find the value of x = X at the point d' arret. This 
may be written as 



-I 



■4k , 

u du 



X=ac 2k 2k + «> 
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where the integral is again taken over the roundabout path of Figure 
2. We may write 

Jr. . u du _„ udu 
r— 4* — t* o — — _ __ 

, M . Ml 4 K = / . U , U I 4 K 

m2 -2, + 2knI 1 + ^ J — w2 -2^ + 2- K >|l + ^ 

where the first integral on the right is taken along the junction line, 
and the second about u = i — \k. The former is approximately 

which is of order /c 2 and vanishes to the present degree of approxima- 
tion. The latter is 



-2« 



u 



, . 1 Kl 
U + l — ^K 

2 u 



1 

1--K 

= - 2m - z N 
l 2i (1 + w) 

"2 K 



. 1 

= — m — - itk. 

The value of X is therefore, 

X = - oe-*" = - o(l - Jitk). 

The point d'arret therefore lies slightly nearer the origin than the point 
of rest. 

To obtain the amount of energy lost by radiation during the swing 
from the starting point of no velocity to the point d'arret we might 
evaluate the integral of K(d?x/dT 2 y as we have evaluated the integrals 
for the time and displacement. But it is simpler to figure the final 
kinetic energy and potential energy and to subtract their sum from the 
initial potential energy |o 2 . The final velocity is v = 4k*, with 
x = X = — a(l — \ttk) . The kinetic energy |« 2 is therefore negligible 
to the order of our approximations. The potential energy is \x 2 = 
§(a 2 — itk) which is less than fa 2 by the amount \-kk. Thus the 
energy radiated during the swing is essentially the same as that 
radiated during one swing on the assumption that the motion is 
practically simple harmonic. 



m 
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We have now concluded a tolerably complete qualitative and 
quantitative discussion of the equation 

dxfePx , 2 \ 2eWdW 

dtW + KX ) = -3?W)' (15) 

which arises by the application of the principle of energy with the 
usual law of instantaneous radiation, and which would arise from the 
integrated form of the energy equation by the same reasoning as is 
applied, after a transformation by integration by parts, namely, that 
if an integral vanishes, its integrand also vanishes, to obtain the 
equation 

dt> +KX J 3<?df {Lb) 

It is probable that if the non-linear equation (15) had been easily 
integrable and had given a motion essentially simple harmonic and 
damped, the transformation to (16) would never have been made. 
And the fact that (15) is difficult of integration, and when actually 
integrated gives an aperiodic discontinuous motion does not seem a 
very satisfactory justification for making the transformation. 

6. Conclusions. One of the first inferences of our work is that, 
whatever we may think of the merits of either of (15) or (16), the 
Abraham-Planck-Rudenberg derivation of (16) is without much mathe- 
matical or physical validity. The derivation of Lorentz and Ritz 11 
is not open to the same objections; for it proceeds directly from the 
retarded potentials to (16) through a series of approximations. The 
law of radiation which this equation states is that the rate of radiation 
is 

dx fd 2 x 2 \ _ 2 e 2 dx d 3 x 
m dt\dJ + KX J ~ 3?dtd¥' 

which for oscillatory motions is nearly proportional to the square of 
the velocity and is physically nearly equivalent to a friction propor- 
tional to the velocity. The average rate of radiation is, of course, 
that which corresponds to that found by Hertz. 

Our work throws some light on the method of presenting the differ- 
ential equation for the oscillator of radiotelegraphy. We may assume 
that the oscillations are damped simple harmonic; we know that the 
simplest mathematical representation of such motion is by an equa- 

11 Lorentz, Theory of Electrons, loc. cit. ; Walther Ritz, Gesammelte Werke, 
401 ff. 
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tion of the second order with constant coefficients, and we may deter- 
mine the logarithmic decrement of this equation so that the average 
loss of energy shall be equal to that determined by the theory of 
Hertz. 12 This is frankly phenomenological and seeks not to go too 
far below the superficial appearances; but it is probably better than to 
delve so deeply that we must sacrifice either our skin or our conscience 
in getting back to the surface. 

We have thus far confined our attention to the oscillator, but simi- 
lar results may be obtained for other systems. For instance, J. J. 
Thomson 13 has proposed as a model of the atom a center repelling a 
corpuscle with a force varying as the inverse cube of the distance (and 
furthermore attracting as the inverse square if the corpuscle lies 
within a tube issuing from the center). To treat the motion under 
the repelling force he writes 

d?x _ ce 
m dP'~ x*' 
integrates the equation and then calculates the radiation by 

COmt J ( W) dL 

The motion is therefore obtained by ignoring the reaction of the radia- 
tion. 

If, however, we write the equation for rate of change of energy at 
each instant, we have an equation of the type 

dr fdfr 1 \ (dhY* 



dt \dt 2 r 3 J """ K Ut 2 



, - v II 
f= ~2T *K<j* + 



-J .4- A I o , 4m> 



For approach to the center, v < 0,/ > 0; but the acceleration becomes 
imaginary when — vr 3 = 4k, and the motion has a point d' arret. 
Motion away from the center does not have this peculiarity. 

The real interest of our work lies, however, not in the critique of 
current "demonstrations," nor in a tour de force integration of a 



12 See, for instance, J. Fleming, Principles of Electric Wave Telegraphy, 
chaps. 3, 5 (1906). 

13 J. J. Thomson, On the Structure of the Atom, Phil. Mag. (6) 26, pp. 792- 
799, Oct. 1913. 
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certain differential equation, nor in pointing out that the integration 
of the equation leads to a physical point d'arret, but in the general 
inferences or alternative inferential possibilities of importance to 
physical science which the discovery of the point d'arret may suggest. 

It is clear that we shall not have to abandon the statement of physi- 
cal problems in terms of differential equations possessing mathemati- 
cally continuous solutions in order to introduce into physics discontinu- 
ities even more startling than those introduced by the adherents of cur- 
rent theories of quanta. The statement of simple cases of rectilinear 
motion in equations of the second degree and second order has led to 
cusps in the space-time locus and, by virtue of the irreversibility of 
time, to physical discontinuities. (Such a conclusion could hardly be 
drawn from Borel's work, for his equation led to no point d'arret.) 

If we regard the existence of the point d'arret as inconceivable, 
we must look about for a method of doing away with it. One possi- 
bility is to accept literally the method of calculation employed by J. J. 
Thomson as cited above, namely, to regard the motion of the electric 
charge as determined by an ordinary equation of mechanics, which 
omits altogether the reaction of the radiation, and to compute the 
radiation subsequently. This would mean that the radiant energy 
did not come out of the field or motion (potential and kinetic energy) 
of the system but was a phenomenon inherently arising when charges 
moved. 14 This would be a shock to those who have implicit confidence 
in the conservation of energy as now generally understood, but would 
not perhaps be so serious to those who take the position of Ritz 15 : 
La localisation de l'energie doit done etre comptee au nombre des con- 
ceptions logiquement inutiles (et peut-etre parfois nuisibles) de la 
theorie . . . Dans le cas le plus general du rayonnement electromagne- 
tique, la conservation de l'energie n'est plus une loi, mais une con- 
vention. 

Another method of escape would be to deny the accuracy of the law 
that the radiation is proportional at each instant to the square of the 
acceleration. The law is indeed only an approximation to the general 
form given by Abraham and Heaviside (with the apparent assumption 
of the quasi-stationary state) and proved by Lewis and me (on the 
basis of the principle of relativity) to be rigorously true for the point 
electron moving with any velocity and acceleration, 16 — for with the 

14 This may remotely suggest Bateman's double-barreled gun in his Corpus- 
cular Radiation, Phil. Mag. (6) 26, pp. 579-585, Oct. 1913. 

15 Gesammelte Werke, pp. 343, 345. 

16 Wilson and Lewis, loc. cit., (5). 
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usual definitions of field and energy the law becomes merely a geometric 
theorem in four dimensional non-Euclidean geometry. It is reason- 
able, however, to suppose that the use of the general form of the law 
could not remove the difficulty; for the point d'arret arises under 
conditions of relatively minimal velocity. Moreover, a qualitative 
discussion of the equation which arises under the general law indicates 17 
that the difficulty is not removed. 

It is tolerably clear that we cannot avoid the point d'arret and main- 
tain simultaneously the point electron, ordinary mechanics, and ordi- 
nary electromagnetic theory. This conclusion may be welcome to the 
school that follows Bohr in abandoning both mechanics and electro- 
magnetic theory 18 as hitherto understood even when transformed as . 
suggested by the principle of relativity. And we certainly do not wish 
to imply that this abandon may not work out satisfactorily. 

The conservative, however, will naturally try to work out of the 
present difficulty by abandoning the point electron, especially as elec- 
trons are generally supposed to have magnitude. It is a bit hard to see 
why the assumption of a very small finite size for the electron must 
fundamentally vitiate the reasoning which leads to the law of the square 
of the acceleration; but we shall not go into this question, because we 
are not acquainted with any derivation of this law which does not in 
some form practically assume the electron to be a point. Moreover, 
the work of Lorentz and Ritz, starting with a distributed electron and 
introducing approximations based on the assumption of small velocities 
and accelerations, establishes a reaction proportional to the rate of 
change of the acceleration and thus leads on mechanical principles to 
a law of radiation (at any rate to a loss of mechanical energy) propor- 
tional to the (scalar) product of the velocity and rate of change of accele- 
ration. And it is interesting to note that this reaction and rate of 
change of energy is actually independent of the size or shape assumed 
for the electron, that is, may be assumed to hold for the point electron. 

Howeverso apologetic one may be in regard to (16), that equation 
seems somewhat less in need of apologies than (15) or its derivation 
from (15) in integrated form. 

Massachusetts Institute of Technology, 
Boston, Mass, May, 1914. 

IV Just what assumptions in non-Newtonian mechanics should be made to 
treat this problem it is perhaps difficult to state, and consequently the question 
is here left as a qualitative judgment of the author. 

18 N. Bohr, Phil. Mag., (6) 26 (1913) a series of three articles On the Constitu- 
tion of Atoms and Molecules, pp. 1, 476, 857. 



